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The low-temperature properties of the spin S = 1/2 ladder with anisotropic ferromagnetic legs are 
studied using the continuum limit bosonization approach. The weak-coupling ground state phase 
diagram of the model is obtained for a wide range of coupling constants and several unconventional 
gapless "spin- liquid" phases are shown to exist for ferromagnetic coupling. The behavior of the 
ladder system in the vicinity of the ferromagnetic instability point is discussed in detail. 
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I. INTRODUCTION 

A theoretical understanding of the properties of quan- 
tum spin ladder systems has attracted a lot of current 
interest for a number of reasons: on the one hand there 
is an increasing number of new magnetic materials with 
a ladder-like structure characterized by rich ground state 
phase diagrams^. On the other hand, spin-ladder models 
pose interesting theoretical problems. For example, since 
a spin-S* chain can be described as a 2S'-leg ladder with 
spin S — 1/2, provided the interchain coupling is appro- 
priately chosen^"** the even- or odd-leg ladder systems 
are an excellent demonstration for Haldane's conjecture^ 
as generalized to S* = 1/2 ladders: the antiferromagnetic 
spin ladder with an even number of legs corresponds to 
a spin chain with integer spin and is predicted to have a 
gap, while a ladder with odd number of legs has a gapless 
excitation spectrum. The two-leg antiferromagnetic lad- 
der is presumably the simplest spin system which allows 
to follow the continuous evolution between spin S* = 1/2 
and S = 1 antiferromagnetic chains nearly exactly^'^. 

The two-leg ladder model has been the subject of con- 
siderable theoretical interest ^"'^^ , most of this work, how- 
ever, concentrated on isotropic or weakly anisotropic an- 
tiferromagnetic chains coupled by an interchain exchange 
of arbitrary sign. Ladder systems with ferromagnetic legs 
are much less studied although, as will be shown here 
they exhibit new interesting aspects. The ferromagnetic 
ladder model in the vicinity of the ferromagnetic instabil- 
ity point was recently studied by Kolezhuk and Mikeska 
within the framework of quasiclassical analysis based on 
a nonlinear cr-model approach^'*. In this work we extend 
these studies addressing the problem of ref. 24 in the 
extreme quantum limit of spins 1/2 and investigate the 
weak-coupling ground state phase diagram of a S* = 1/2 
ladder system with ferromagnetic legs and anisotropic in- 
terleg exchange using the continuum limit bosonization 
approach. 

The Hamiltonian of the model under consideration is 
given by 

H^Hl^ + Hl^+H^, (1) 
where the Hamiltonian for leg a is 



N 
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^leg ~ 'J Waj'-'ctj + l ^ a,j a,j + l 

J = l 

+ ASl^Sl^+,), (2) 
and the interleg coupling is given by 

N 
N 

+ JIJ2SI^SI^. (3) 

i=i 

Here S^''^'^ are spin 5=1/2 operators at the j-th 
rung, the index a = 1, 2 denotes the ladder legs. The 
intraleg coupling constant is ferromagnetic, J > 0, and 
therefore the limiting case of isotropic ferromagnetic legs 
corresponds to A = 1, while the case of isotropic antifer- 
romagnetic legs is obtained for A = — 1. Since we study 
the ground state phase diagram of the ladder system with 
ferromagnetic legs we will restrict ourselves to consider 
A > 0, including | A| ^ 1 in the limit of legs with strong 
in-plane anisotropy. 

The outline of the paper is as follows: In section II we 
derive the bosonized formulation of the model in the con- 
tinuum limit. In section III wc discuss the weak coupling 
phase diagrams of the model for three different cases of 
anisotropic interleg coupling. Finally, we conclude and 
summarize our results in section IV. In the appendix we 
present the spin-wave approach to study the transition 
line related to the ferromagnetic instability in the sys- 
tem. 



II. BOSONIZATION 

In this section we derive the low-energy effective field 
theory of the lattice model Eq. (1). Since the weak- 
coupling bosonization approach to the ladder models is 
based on the perturbative treatment of the interchain 
couplings^^''* we start from the bosonization of separate 
spin S — 1/2 ferromagnetic XXZ chains. 
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A. Separate chains 

The anisotropic spin 5=1/2 Heisenberg chain with 

|A| < 1 is known to be critical. The long wavelength ex- 
citations are described by the standard gaussian theory"^® 
with hamiltonian 



0? 



(4) 



4'{x) and 9{x) are dual bosonic fields, dtcj) = ud^O, and 
satisiy the following commutational relation 
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Note that the first and the second terms in Eqs. 
(9), (10) are hermitian because of (5). Furthermore : ... : 
denotes the normal ordering with respect to free bose 
system (4), a is the leg index. The non- universal real 
constants a, h and c depend smoothly on the parame- 
ter A, are of the order of unity at A = 0^°'^^ and are 
expected to be nonzero at arbitrary A < 1. 



u stands for the velocity of spin excitation and is fixed 
from the Bethe ansatz solution as'^^ 



B. Coupled Spin-1/2 Chains 



K 



(6) 



where K is Luttinger liquid parameter known from com- 
parison with the exact solution of the XX Z chain: 



K = 



2 arccos A 



(7) 



Thus the parameter K increases monotonically along 
the XX Z critical line — 1 < A < 1 from its minimal 
value i4'=l/2atA = — 1 (isotropic antiferromagnctic 
chain), is equal to unity at A = (the XY chain) and 
ii'^ooatA— >1. AtA = l the spin excitation velocity 
vanishes, u = 0. This corresponds to the ferromagnetic 
instability point of a single chain. 

To obtain the bosonized version of the ladder Hamil- 
tonian we need the explicit bosonized expressions of 
the spin operators. The bosonization procedure for the 
spin S = 1/2 Heisenberg chain is reviewed in many 
places'^ However, since we consider the ladder model 
with ferromagnetic legs, our bosonization conventions re- 
quire some comments. The unitary transformation 



a,] 



i-iys; 



a,3 



(8) 



changes the sign of the intrachain transverse exchange 
and maps the Hamiltonian (1) to the Hamiltonian with 
antiferromagnctic legs. To maintain the ferromagnetic 
character of the in-plane correlations in the bosonization, 
it is convenient to implement the multiplicative factor 
(— 1)-', introduced by the unitary transformation (8), 
directly in the bosonized expressions for the transverse 
components of the spin operators. Using the standard 
bosonization formulas^^"^^ we obtain: 
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For coupled 5=1/2 chains we start with two critical 
S = 1/2 Heisenberg chains and treat the inter leg coupling 
as a perturbation assuming | Jj^ | , | \ ^ J. Therefore 
we start with two free bose field Hamiltonians (4) and 
simply attach a leg index a = 1, 2 to the fields. 

Then we introduce the symmetric and antisymmetric 
combinations of the bosonic fields: 



9± = 



2A 



(12) 



where 



and 



A± 



Jeff = J 
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Using Eqs. (9)-(ll), we finally obtain the following 
bosonic Hamiltonian density: 
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Here 
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u± = - — ~ w 1 ± 

K± = K ■ A±c^ K (It ^ 



(17) 
(18) 



27r Jeff , 

and we have introduced the following coupling constants: 
Jl = Jlh, (19) 
for A = and otherwise 



1^ 



(20) 
(21) 



with some positive constants of proportionahty which 
cannot be fixed by symmetry arguments (contrary to the 
constant appearing in Eq. (19)). 

In deriving (13), a term ^ cos \J j^O- cos y^SnK^cj)- 

which is strongly irrelevant at A > (ferromagnetic legs) 
was omitted. Thus our approach is tailored to cover fer- 
romagnetic intraleg coupling and can be applied to anti- 
ferromagnetic intraleg coupling only for |A| <C 1. 



C. The effective continuum-limit model 

At J]_ = = the Hamiltonian (13) describes two 
independent Gaussian fields, i.e. two gapless fields, each 
describing a critical spin S=l/2 Heisenberg chain. Let us 
first address the question whether the interleg exchange 
leads to the dynamical generation of a gap in the exci- 
tation spectrum. In the case of anisotropic interleg ex- 
change it is rather instructive to study the effect of the 
longitudinal ( J^) and transverse (J^^) part of the inter- 
leg coupling separately. 

At — the effective theory of the original ladder 
model is given by two decoupled quantum sine-Gordon 
models 



Tie// = + Ti 



(22) 



where 



7i = u± I dx 



+ — — cos \JS'!TK±(f)± {x) 



(23) 



The two SG models respectively describe the symmetric 
((/>+) and antisymmetric (0_) degrees of freedom. 

The bare values of the dimensionless coupling con- 
stants M± and K± are known only in the weak-coupling 
limit I J^l/J, |A| <C 1 where they have the values: 



nJ 
2A 

i^± = l + — 



(24) 

(25) 



The scaling dimensions of the cosine terms in (23) 
2K± ~ 2. Therefore, in the weak-coupling 



limit, both SG models have marginal dimension and de- 
tails of their behavior should be determined within the 
framework of the renormalization-group analysis. How- 
ever, a rather straightforward estimate indicates that at 
A ^ 4J^/J we have K± > 1 and the cosine terms are 
irrelevant. Therefore in this case one concludes that the 
effective model reduces to the theory of two independent 
Gaussian fields ^± . The effect of the interleg coupling is 
extremely weak and is completely absorbed in the renor- 
malized values of the spin-liquid parameters K± char- 
acterizing respectively gapless symmetric and antisym- 
metric modes. However, it is important to note, that 
in the vicinity of the single chain ferromagnetic instabil- 
ity point, at A 1, the effective bandwidth collapses, 
Jeff ~ 2(1 — A) — » 0. Therefore in this limit the effect of 
the interleg coupling becomes very strong. This implies 
subtle effects to be discussed later. 

The transverse interleg exchange (J^^) leads to 
the appearance of the strongly relevant operator 

cos with the scaling dimension ct_y = 

(2K^)^^ < 1/2 in the theory. Therefore, the antisym- 
metric sector is gapped at arbitrary J^ ^ 0. Fluctua- 
tions of the field 9- [x) are completely suppressed in this 
sector and 0^{x) is condensed in one of its vacua. The 
vacuum expectation value of the cosine term is 



(cosy2^riF^6l_) =7 



(26) 



with 7 ~ {\J±^\I Jeff) ^~ ~ ^ 1 in weak coupling and is 
of the order of unity at \J^\ > J- 

Therefore, the condensation of the field 9^ strongly 
influences the coupling between the symmetric and anti- 
symmetric modes induces by Hf^^. Taking into account 
that the fluctuations of the fleld 9- are stopped, one eas- 
ily flnds that at J^ ^ infrared behavior of the sym- 
metric field is governed by the following "effective" sine- 
Gordon theory 

Htff = u+jdx \^-[{dJ+f + {dMx)f] 



where 



+ cos J^TTK^(i)+ ix) 
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(27) 
(28) 



are 



This mapping of the initial spin 5=1/2 ladder model 

onto the quantum theory of Bose fields described in terms 
of an "effective" sine-Gordon (SG) Hamiltonian (23) or 
(27) will allow to extract the ground state properties of 
the 5=1/2 ladder using the far-infrared properties of 
the quantum SG theory. 



D. The RG analysis 

The infrared behavior of the SG Hamiltonian is de- 
scribed by the corresponding pair of renormalization 
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group (RG) equations for the effective coupling constants 
/C(0 and M{1) 



dM{l) 

dl 
dK{l) 

dl 



= -2{]C{l)-l)M{l) 



(29) 



where I = In(ao) and the bare values of the coupling con- 
stants are /C(Z = 0) = K a.nd M{1 ^ 0) = M. The pair 
of RG equations (29) describes the Kostcrlitz-Thouless 
transition'*^. The flow lines lie on the hyperbola 



(30) 



and exhibit two different regimes depending on the rela- 
tion between the bare coupling constants -^'^ (see Fig.l): 



MA 




2(K-1) 



FIG. 1. Renormalization-group flow diagram; the arrows 
denote the direction of flow with increasing length scale. 

Weak coupling regime. For 2{K — 1) > \M\ we are 
in the weak coupling regime: the effective mass A4 ^ 0. 
The low energy (large distance) behavior of the corre- 
sponding gapless mode is described by a free scalar field. 

The vacuum averages of exponentials of the corre- 
sponding fields show a power-law decay at large distances 



iK<l>i_0) „-iK4>{r) 



iKe{0) -iKe{r) 



(31) 



where K* is the fixed-point value of the parameter K 
determined from the Eq. (30). 

Strong coupling regime. For 2{K — 1) < \M\ the sys- 
tem scales to strong coupling: depending on the sign of 
the bare mass M , the renormalized mass M is driven to 
±oo, signaling a crossover to one of two strong coupling 
regimes with a dynamical generation of a commensura- 
bility gap in the excitation spectrum. The flow of | | 
to large values indicates that the AicosV^TT K(l> term in 
the sine-Gordon model dominates the long-distance prop- 
erties of the system. Depending on the sign of the mass 
term, the field (j) gets ordered with the expectation values 



^/tt/SK at M > 







at M < 



(32) 



Using this analysis for the excitation spectrum of the SG 
model and the behavior of the corresponding fields, Eqs. 



(31, 32), we will now discuss the weak-coupling phase 
diagram of the spin S = 1/2 ferromagnetic ladder model 
Eq. (1). 



III. PHASE DIAGRAMS 

In this section we discuss separately the ground state 
phase diagram of the ferromagnetic ladder coupled only 
by the longitudinal part of the interleg spin exchange 
(subsection A) , coupled only by the transverse part of the 
interleg spin exchange (subsection B) and by an isotropic 
interleg coupling (subsection C). At this point we note 
that from the structure of the interaction Hamiltonian 
Eq. (3) follows that the phase diagrams for case (A) 
and case (B) will be symmetric with respect to the lines 
Jj_ = since a change of sign in Jj_ leads to a unitary 
equivalent Hamiltonian. This is in contrast to case (C) 
where this unitary equivalence does not exist. 



A. Chains coupled by the longitudinal part of the 
interleg exchange 

In this subsection we consider the weak-coupling phase 
diagram of the spin S=l/2 ferromagnetic ladder model 
(1) coupled by a weak longitudinal intrachain exchange 
{J^ = 0, Jf ^0). The bosonized version of the model 
Heff = H'^+Tl~ where Tt^ are given by Eq. (23) and the 
bare values of the corresponding dimensionless coupling 
constants are given by (24)- (25). By inspection of the 
initial values of the coupling constants one easily finds 
that: 

• At A < both the symmetric and the antisymmet- 
ric sectors are gapped (except for =0). 

• At A > the symmetric sector is gapped for 
J^/J > 2 A > while the antisymmetric sector 

is gapped for J^/J < — 2A < 0. 

This determines the following three distinct sectors of 
the phase diagram as traced already in the RG analysis 
(see also Fig. 2 below): 

• Sector A: A < corresponds to the phase with 

gapped excitation spectrum; 

• Sector B: A > and \J]_\ > 2JA corresponds to 
the phase characterized by the one gapless and one 
gapped mode in the excitation spectrum. In par- 
ticular at > the symmetric mode is gapped 
whereas the antisymmetric mode is gapless and 
vice- versa at Jj_ < 0; 

• Sector C: A > and \J^\ < 2JA corresponds to 

the phase where both modes are gapless. 

As we show below, the same phases are present in the 
strong coupling regime. The only phase which is missed 
in the weak-coupling RG analysis is the ferromagnetic 
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phase; it appears only in the strong coupling regime at 
A ~ 1 or at A < 1 but I J^l ~ 1/A » 1. 

To clarify the symmetry properties of the ground states 
of the system in the different sectors we study the large- 
distance behavior of the longitudinal 



and the transverse 



(33) 



(34) 
P) and 



spin-spin correlation functions for intraleg (a 
interlcg (a ^ [3) spin pairs. 

Using the results for the excitation spectrum and the 
behavior of the corresponding fields in the gapless and 
gapped phases, Eqs. (31)-(32), and the expressions for the 
corresponding correlation functions from bosonization, 
we now discuss the characteristics of the various phases 
in the different sectors of the weak- coupling ground state 
phase diagram. 

In the sector A ( < 0) the vacuum expectation val- 
ues of the fields are: = ^^tt/8K^ and = 0. 
Ordering of the 0± suppresses transverse spin correla- 
tions, while the longitudinal correlations are given by 



Km 



{—ly ■ constant. 



Therefore at A < and < 0, the long-range ordered 

(LRO) antifcrromagnctic phase with inphasc spin order- 
ing on the rungs is realized in the ground state of the 
system. 

In the Sector Al (J^ < 0) the vacuum expectation 
values of the fields are given by (0+) = and (</)_) = 
^/^^/8KZ. This immediately implies that in this sector 

K^pir) - (-1)"+^'^ . (-l)-^ . constant. 

Therefore, at A < and > the LRO antiferro- 
magnetic phase with antiphase intrarung spin ordering is 

realized in the ground state of the system. 

In the sector B (Bl) the antisymmetric (symmet- 
ric) field is gapped with the vacuum expectation value 
((/)_) = ({4>+) = 0). However, as can be seen from weak 
coupling RG analysis as well as from the strong coupling 
effective spin 1/2 model (see below), at ^ the line 
A = is the phase transition line along which the gapped 
at A < symmetric (antisymmetric) mode becomes gap- 
less. Therefore in the sector B (Bl) the gapless degrees 
of freedom corresponding to the symmetric (antisymmet- 
ric) mode are described by the free Bose field system with 
the fixed-point value of the parameters K^. Using Eq. 
(30) and the bare values of coupling constants (24)-(25) 
it is straightforward to show that at |A|, | J^/J| <C 1 



1 



1 + -J2A(2At JI/J). 



Note that at A = the fixed-point values of the spin- 
liquid parameters are = 1 while at = (see Eq. 
(18)) Kj_ = K. Therefore we conclude that along the 



line A = the gapless sector in the system is identical to 
a single isotropic spin 5 = 1/2 Heisenberg chain, while 
along the line J\_ = we reach the limit of two decoupled 
spin S = 1/2 Heisenberg chains. 

The very existence of a gapped excitation mode ac- 
companied with ordering of the field (p- (or (/)+) implies 
suppression of the transverse correlations. On the other 
hand the presence of the gapless excitation mode leads to 
the power law decay of the longitudinal spin correlations. 
Therefore we obtain, that in the sector B 



while in the sector Bl 



27rr2 



-I- 



(-1)^ 



We denote this phase as the spin liquid I phase. It is 
interesting to note that in sector B the following operator 
shows quasi long range behavior: 



c.4(5+(r)5+(r))5r(0)52-(0)) 



(35) 



{S" + 5*2 )^ for the S = ^ ladder corresponds to the op- 
erator (S^y in the 5 = 1 chain and we therefore identify 
sector B with the XY2 phase for the S = 1 chain as 
described in ref.[ 2]. 

With increasing interleg ferromagnetic coupling we 
reach the line A = \J^\/2J which marks the transition 
into the phase where both fields are gapless. In the sec- 
tor C of the phase diagram the system shows proper- 
ties of two almost independent spin 5=1/2 anisotropic 
Heisenberg chains with dominating ferromagnetic cou- 
pling. The transverse correlations in this phase are given 
by 



(_!)'■ . ^-{K'++K'_+l/4Kl+l/4K'_ ) 



(36) 



where 6a, i3 is the Kronecker symbol. The longitudinal 
correlations decay faster. In particular the intraleg lon- 
gitudinal correlations are given by 



aa 



(r) 



A': 



A: 



27rr2 



+ {-iy-r 



(37) 



The transverse interleg correlations are strongly sup- 
pressed in this phase, while the longitudinal part of the 
interleg spin-spin correlations is given by 



K* 



27rr2 



(38) 



This phase we denote as the spin liquid II phase. 

Although the analysis as considered above is formally 
valid in the weak-coupling limit (A, | <C J) we can 
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estimate the upper boundary for the spin liquid //phase, 
in the vicinity of the single chain ferromagnetic instabihty 
regime using the dimensionahty analysis. We determine 
the instability curve corresponding to the transition into 
the gapped phase from the condition K± = 1, where 
the scaling dimension of the corresponding cosine term 
d± = 2K± = 2. After some simple algebra one easily 
obtains that at > the field (p- is gapless, while the 



field becomes massive for 

Jl> Jl = 2-KU- 



K-1 
1^ 



(39) 



For A = 1 — e with e <C 1 , which implies 1/K • 
equation (39) takes the following form 



2e<Cl 



and the pseudospin operators are 



2 ("^n,! 



2 > 



2 ' 



^n,2) ■ 



In agreement with the weak-coupling bosonization 
analysis, at A = (XY-legs) the system is equivalent 
to the S* = 1/2 isotropic antifcrromagnctic chain. For 
arbitrary A < {)^lff > ^eff) ^P^'^ chain given by 
the Hamiltonian (42) is in the gapped Neel phase. This 
phase corresponds to the LRO AFM intcrleg ordering 
with interleg phase shift equal to zero. At 

< A < J/\Jl\ 



J^(A) =4Je 1- 



2V2e 



(40) 



Therefore in the vicinity of the single chain ferromag- 
netic instability point, at 1 — A <C 1, the spin liquid I 
phase with only one gapless (here antisymmetric) mode 
reenters the phase diagram at Jj_ > J^(A). (We note 
that the amplitude of the cosine term in the limit of the 
single chain ferromagnetic instability point is not deter- 
mined exactly, so the phase transition line determined by 
the dimensional analysis is of qualitative nature in this 
limit). For < the analysis is done in exactly the 
same manner with symmetric and antisymmetric modes 
changing roles. 

At = and A > 1 each of the decoupled legs 
is unstable towards the transition into a ferromagnetic 
phase. At 7^ 0, we can address the problem of the 
ferromagnetic instability in the ladder system studying 
the velocity renormalization of the corresponding gapless 
excitations. In analogy with the single chain case we 
mark the transition into the ferromagnetically ordered 
phase at u± = 0. Using Eqs. (6) one finds that the 
ferromagnetic transition takes place at 



J^M(A)=4Je. 



At I I 3> J the boundary of the ferromagnetic insta- 
bility can be established from the large rung coupling ex- 
pansion approach. Let us first consider the case of strong 
ferromagnetic intrarung interaction < 0. In this limit 
a large gap of order | Jj^| exists in the one-magnon ex- 
citation spectrum. Projecting the system on the sub- 
space excluding antiparallcl orientation of spins within a 
given rung, in the second-order perturbation expansion 
withr respect to J'^/\Jj_\ and up to the additive constant 
Eo = —Nq\J^\ wc obtain the following effective spin-1/2 
XXZ spin chain Hamiltonian 



(42) 



< Af f , ) the spin chain (42) is in a gapless 



planar XY phase, corresponding to the "spin liquid I" 
phase of the bosonization studies and finally at 



A> J/IJll 

i^lff < —Xgff) the transition into the completely polar- 
ized ferromagnetic phase takes place. 

In the case of strong antiferromagnetic interleg cou- 
pling J]^ ^ J > analysis is similar. In this case the 
intrarung ordering of spins is antiferromagnetic. Project- 
ing the system on the subspace excluding parallel orien- 
tation of spins within the same rung, and introducing a 
new set of spin operators 



,2 ' 



2 ' 



in the second-order with respect to J'^/Jj_ we once again 
map the initial ladder model onto the theory of an 
anisotropic spin 1/2 Heisenbcrg chain (42). One can per- 
form the analysis as discussed above, however the ferro- 
magnetic ordering in terms of the effective 5=1/2 chain, 
(41) at Jj^ > corresponds to an intcrleg ferromagnetic order- 



ing with a phase shift of tt of the order parameter along 
the rung. 



FERROMAGNET 
(Inphase interchain order) ^ 




FERROMAGNET 
^ (Antipliase intercliain order) 



B(XY2) 

Spin Liquid I 



Liquid 

II Jjy---- ------ TBI" 

Spin Liquid I -- 



LRO AFM A 
(Inphase interchain order) 



Al LRO AFM j//J 
(Antiphase interchain order) 



where 



\J1\' 



Kff - 



2JA 



(43) 



FIG. 2. The groundstatc phase diagram of the two-leg lad- 
der with a longitudinal JxS|,i>Sj_2 coupling between legs. For 
details see text. 
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The results obtained within the bosonization approach 
together with the results from the strong coupling expan- 
sion allow to draw the following phase diagram of the lad- 
der with a longitudinal inter leg coupling (see Fig. 2). 
At A < the phase diagram consists of two gapped 
phases describing respectively long range ordered Neel 
antiferromagnetic phases with gapped excitation spec- 
trum and inphase (at Jj^ < 0) or (antiphase at J]_ > 0) 
ordering of spins within the same rung. The line A = 
marks the transition into the Spin Liquid I - phase char- 
acterized by a gapless excitation spectrum and power law 
decay of the spin-spin correlation functions. The criti- 
cal indices for the decay of the corresponding spin-spin 
correlations in the Spin Liquid I phase are ji ~ 1. In 
the case of strong interleg exchange \Jj_\ » J, further 
increase of the interleg ferromagnetic exchange A leads 
to the transition at Ac ~ J/\J]_\ into the phase with 
ferromagnetically ordered legs. However, in the weak- 
coupling case, at | Jj^ | ^ J, an increase of the parameter 
A at given Jj^ leads to the transition into the Spin Liquid 
//at Ac(i) = I 1/2 J. The Spin Liquid II phase is char- 
acterized by a gapless excitation spectrum and power-law 
decay of the spin-spin correlation functions with critical 
indices 7^ ~ 2. This transition marks the development of 
the regime dominated by intraleg coupling, whereas the 
interleg longitudinal exchange plays only a rather mod- 
erate role. However, with further increase of the intraleg 
ferromagnetic exchange, in the vicinity of the ferromag- 
netic instability line the Spin Liquid II phase becomes 
unstable and the system reentres into the Spin Liquid I 
phase. This reentrance effect is connected with a sharp 
reduction of the bandwidth in the vicinity of the fer- 
romagnetic transition and a subsequent increase of the 
potential energy of the interleg coupling. Therefore, just 
before the transition into the ferromagnetically ordered 
phase, the short range interleg fluctuations get stopped, 
and as in the case of the strong inrtarung coupling, the 
Spin Liquid I - phase is unstable toward the transition 
into the phase with ferromagnetically ordered legs. 

However, since the transition into the ferromagnetic 
phase is a typical finite bandwidth effect, the parame- 
ters determined quantitatively within the bosonization 
(i.e. infinite band) approach strongly depend on the 
way of regularization of the continuum theory on small 
distances. Therefore, it is useful to determine the low- 
est boundary of the ferromagnetic phase on the phase 
diagram, starting from the ferromagnetically ordered 
phase and using the standard spin- wave analysis (see ap- 
pendix). At \ JW ^ J this approach gives Jp^^ — 2(/e. 
This discrepancy clearly is the result of the linearized 
expressions used for the parameters K±, Eq. (18). How- 
ever, as long as the multiplicative renormalization used 
here for the parameters u± and K±, Eqs. (17)-(18) re- 
mains valid, the scenario discussed above, where the spin 
liquid I phase is unstable towards the transition into the 
phase with ferromagnetically ordered legs, remains qual- 
itatively plausible. For a more quantitative description, 
sufficiently detailed numerical studies of this sector of the 
phase will be very helpful. 



To conclude this subsection we note that the groimd 
state phase diagram of the ferromagnetic ladder system 
coupled only by the longitudinal part of the spin-spin ex- 
change interaction exhibits a rather rich phase diagram 
which consist of LRO AFM phases, a spin liquid phase 
with one gapped and one gapless mode, a spin liquid 
phase with two gapless modes and a phase with ferro- 
magnetically ordered legs. 



B. Chains coupled by the transverse part of the 
ladder exchange 

In this subsection we consider the case of two critical 
Heiscnberg chains coupled by a transverse interleg ex- 
change interaction = and ^ 0. The particular 
aspects of this limiting case are the following ones: 

• The antisymmetric mode is gapped at arbitrary 
J7 ^ 0. 

• The low energy properties of the system are deter- 
mined by the behavior of the symmetric field. 

• The infrared properties of the symmetric field are 
determined by the subtle coupling between the 
symmetric and antisymmetric modes. 

We start our analysis from the limiting case of weakly 
anisotropic XY chains, coupled by the weak interleg 
transverse exchange, assuming \A\,\J^\/ J <i; 1. At 
7^ the antisymmetric mode is gapped and the dual 
antisymmetric field is " pinned" with vacuum expectation 
value 



{0-) = 



^J-kK_/2 at J7 > 







at J7 < 



(44) 



Behavior of the symmetric field is governed by the SG 
Hamiltonian (27). The standard RG analysis gives that 
the symmetric mode is gapped at 



A < Aci 



1_ 
4J 



\J+- 



(45) 



Therefore at A < Ad the excitation spectrum of the 
system is gapped. The dynamical generation of a gap in 
the symmetric mode leads to condensation of the field 
0+ with a vacuum expectation value (</>+) = 0. Since the 
dual component of the antisymmetric field is "pinned" 
with vacuum expectation value given by (44), the so- 
called "disordered" phase-^ is realized in the ground state. 
At J^^ > 0, spins on the same rung form a singlet and 
the ground state corresponds to the state with a singlet 
pair on each rung. There is no correlation between spins 
along the ladder. In the case of ferromagnetic coupling, 
at J^^ < spins on the same rung form a state cor- 
responding to the 5^ = component of the triplet (an 
"asymmetric triplet" pair) and the ground state corre- 
sponds to the state with an "asymmetric triplet" pair on 
each rung. In analogy to the phases of the 5=1 chain 
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as discused in^ we denote this phase as "anisotropic large 
D phase". 



Anisotropic 
Xarge D' phase 



Singlets 




Anisotropic 
Xarge D' pliase 



Rung J^'/J 
Singlets 



FIG. 3. The groundstate phase diagram of the two-leg lad- 
der with transverse coupling between legs. 

For A > Aci the system is in the phase where the sym- 
metric mode is gapless. Since the antisymmetric mode 

is gapped the alternating part of the spin-spin correla- 
tions is exponentially small, while the smooth part shows 
a power law decay at large distances. In particular, the 
in-plane correlations are given by 



K^^{r) ~ constant ■ I — 



\JT\ 



a+l3 



-1/4KI 



while the longitudinal correlations decay faster 



2^ 



(46) 



(47) 



As follows from Eq. (46) the Une J^^ = marks the 
transition from a regime with ferromagnetic interleg or- 
der into the regime with antifcrromagnctic intcrlcg order. 
In the vicinity of this critical line = the gap in the 
antisymmetric mode (M_) is tiny, therefore the correla- 
tion fimctions given by Eqs. (46)- (47) arc valid only for 
distances r ^ Lc ^ 1/M_. However, at distances r <^ 
fluctuations of the antisymmetric mode are strong, and 
the behavior of correlation functions is the same as in the 
spin liquid II phase Eqs. (36)- (38). Following Schulz^ 
who has discussed a similar phase in the context of the 
spin S = 1 chain we denote this phase as a spin liquid 
XYl phase. 

Let us now discuss the phase diagram of the model in 
the vicinity of the single chain ferromagnetic instability 
point A = 1. As A — > 1, the effective coupling constant 
behaves as 



27rM 



2J 



where e = 1 — A. A rough estimate of the renormalization 
of the velocity of the symmetric mode excitations u+ and 



of the spin-liquid parameter 
order gives 



at A ~ 1, in second 



u+=u(l + XK{J+_/jf^ , 



(48) 
(49) 



where A is a nonuniversal constant of the order of unity. 
As follows from Eqs. (48/49) the strong effective trans- 
verse coupling reduces the tendency towards ferromag- 
netic ordering and leads to the transition into the gapped 
phase at 

J+- > J(l-A)i. 
Equivalently we have 

To summarize this subsection we note that the ground 
state phase diagram of the ferromagnetic ladder system 
coupled by the transverse part of the spin-spin exchange 
interaction only also exhibits a rich phase diagram which 
consist of the " disordered rung-singlet" and " anisotropic 
large D" phases, the easy-plane gapless XYl phase and 
the "stripe" ferromagnetic phases with dominating intra- 
leg ferromagnetic ordering. 



C. Chains coupled by isotropic interleg exchange 

In this subsection we consider the weak-coupling 

ground state phase diagram of the model (1) in the case 
of an isotropic interleg exchange = =■ J±- 

In this case the behavior of the antisymmetric sector 
is completely similar to the above considered case of the 
ladder with transverse exchange: the antisymmetric field 
is gapped and the vacuum expectation value of the dual 
field 6- , depends on the sign of exchange and is given by 
Eq. (44) after the substitution J^^ J±. 

The far-infrared properties of the symmetric field are 
governed by the effective SG Hamiltonian (27) with the 
bare values of the model parameters given by and 



M+, = 



-J± (l + S) , 



(50) 



(where 5 is a, nonuniversal positive number) 

The resulting asymmetry of the model is clearly seen: 

• at J7j_ > 0, the antifcrromagnctic interleg exchange 
reduces K+ and increases M^^j: and therefore sup- 
ports the tendencies towards development of a gap 
in the excitation spectrum 

• at < 0, the ferromagnetic interleg exchange in- 
creases while with increasing |j7j_| the param- 
eter M^j^j: ~ J7j_ (1 — (5) ^ ; therefore we expect 
an enlargement of the gapless section in this case. 
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Wc start our analysis from the limiting case of weakly 
anisotropic XY chains assuming |A|,|j7^|/J ^ 1. At 
A = we have K = 1 and the system shows a gap in the 
excitation spectrum at J7j_ > and is gapless in the case 
of ferromagnetic inter leg exchange J'± < 0. Therefore at 
A = 0, with increasing ferromagnetic inter leg exchange 
{J± < 0, \J'±\ oo) the system continuously evolves 
into the limit of the 5 = 1 XY model, which is known to 
be gaplcss'*^^'*^. In the case of antiferromagnetic interleg 
exchange J7j_ > the symmetric mode is unstable to- 
wards the Kosterlitz-Thouless type transition associated 
with the dynamical generation of a gap in the excitation 
spectrum. The weak-coupling RG analysis tells, that at 
A ^ and J71 > the gapless XYl phase is realized for 



A> A 



cl 



J± 

2J 



(l + <5) 



(51) 



whereas in the case of ferromagnetic interrung exchange 
i7_L < it is realized for 



A>A;i 



6J± 
' 2J 



Therefore, from the RG studies we obtain that the gap- 
less XYl phase is stable in the case of ferromagnetic ex- 
change. At J'± > it is unstable towards the transition 
into the gapped rung-singlet phase. At JTI < the gap- 
less XYl phase penetrates into the A < sector of the 
phase diagram. However, since M^ff with increasing 
ferromagnetic exchange, at |j7i| » J the gapless phase 
on the antiferromagnetic side (A < 0) of the phase dia- 
gram shrinks up to a narrow stripe which exponentially 
disappears as |J71|/J — > oo. 

With A ^ 1 the gapless XYl phase becomes unsta- 
ble towards transition into the ferromagnetically ordered 
state. Following the route developed before, we find that 
at A = 1 — e and antiferromagnetic interleg exchange, 
J7j_ > 0, the reentrance of the gapped rung-singlet phase 
takes place at 



Ae2 = 1 - 



J± 
4J 



+ 



J± 
4J 



3/2 



(52) 



Thus, in agreement with the quasiclassical studies^^, 
wc obtain that two almost ferromagnetically ordered 
chains coupled by an isotropic interleg exchange are 
unstable towards formation of the gapped rung-singlet 
phase at J± > J1 > 0, where — > as A ^ 1. How- 
ever, in contrast to the quasiclassical case, J'^^ increases 
linearly with e in the quantum spin-ladder case. 

In the case of ferromagnetic interleg exchange, J'± < 0, 
the gapless XYl phase becomes unstable towards the 
transition into the ferromagnetically ordered phase when 
A increases towards 1. In this case the spin- wave ap- 
proach (see appendix) gives that the boundary between 
the XYl and the ferromagnetic phase is A = 1. 

We summarize our results considering the phase di- 
agram of the ladder with ferromagnetic legs and an 
isotropic interleg exchange in Fig. 4. 



FERROMAGNET 




Anisotropic 
Haldane pliase 



Rung 
Singlets 



FIG. 4. The ground state phase diagram of the two-leg 
ladder with an isotropic interleg coupling. 



IV. CONCLUSIONS 

Wc; have studied the ground state phase diagram of 
the 5 = 1/2 ladder with ferromagnetically interacting 
legs using the continuum limit bosonization approach. 
The phase diagrams for the extreme anisotropic inter- 
chain coupling cases (Ising and XY interleg exchange) 
as well as for the SU{2) symmetric case were obtained. 
These phase diagrams exhibit a number of interesting 
phases, gapped as well as gapless; some of these are fa- 
miliar from well known ID models (rung singlet phase, 
anisotropic Haldane phase, ferromagnetic and large D 
phase) , in addition we describe here for the first time less 
conventional phases for ladders: the spin liquid phases 
with (i) one gapless and one gapped mode (including 
the known XYl and XY2 phases) and (ii) two gapless 
modes. We have shown moreover that the gapped rung 
singlet phase found semiclassically to appear for an arbi- 
trarily small isotropic antiferromagnetic interaction be- 
tween ferromagnetic legs^** continues to exist for 5=1/2 
ladders and xy— like interactions and actually extends to 
small values of A. 

The neighborhood of the single chain ferromagnetic 
instability point turned out to be of particular inter- 
est. We investigated the behavior of the system in this 
regime using the multiplicative regularization scheme. 
This scheme allows to extend the bosonization formal- 
ism to the limit when the bandwidth of the single chain 
excitations collapses and leads to the result that upon 
increasing the strength of ferromagnetism A at any mod- 
erate fixed interleg interaction a sequence of two phase 
transitions occurs before the system enters the final fer- 
romagnetically ordered phase. 

Preliminary investigations show that the system con- 
sidered here displays additional interesting aspects when 
an external magnetic field (both longitudinal and trans- 
verse) is applied. These investigations will be reported 
in a subsequent publication. 
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A = l + ^(JI-Jr). (A3) 

In the particular limit of noninteracting chains ( Jj^ = 
= 0) as well as in the limiting case of the rotation- 
ally invariant interleg coupling (J^ = J^) the critical 
line corresponding to the instability of the ferromagnetic 
phase is given by the condition 

A = 1 . (A4) 



APPENDIX A: FERROMAGNETIC 
INSTABILITY 



In this appendix we consider ferromagnetic interleg 
coupling, assuming J^^, Jj_ < 0, and use the spin wave 
approach to determine the critical line corresponding to 
the ferromagnetic instability in our system. For this pur- 
pose we start from the region of the phase diagram where 
we can safely assume that the ground state is the fully po- 
larized ferromagnetic state (that is A ^ 1 and Jj_ < 0). 
We identify the transition line from the fully polarized 
ground state to some other state as the line of instability 
in the spin wave excitation spectrum. 

Let us denote the eigenstate of the Hamiltonian (1) 
corresponding to the fully polarized (along the Z axis) 
ferromagnetic state by |0). Then 

Sj'^ |0) = 0, for arbitrary j and a. 

It is straightforward to obtain that -ff |0) = Eo\0) where 

Eo = ~{\Jl\+2JA). 

To construct the lowest excitations in the ferromag- 
netic phase we act on the ground state configuration by 
the spin lowering operator S~ ^. Let us denote by 
(|2)^) the state obtained by action of the spin lowering 
operator on the n*^ site of leg 1 (leg 2): 

|l)„^5-i|0), |2)„^5-2|0). 

It is straightforward to solve the coupled system of 

equations of motion in the subspace S^^^ = N — 1 and to 
obtain the following two sets of excitation frequencies: 

a;-(g) = -Jcosg + JA-^(J1- J^^) (Al) 
u;+{q) = - Jcosg + JA - ^ (J^ + J^) . (A2) 

For ferromagnetic interleg exchange {J]_,J'^ < 0) we 
have 

tjJ~{q) < oj+{q) 

and from the instability condition a;~(g = 0) = we 
obtain 
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